Contact angles on heterogeneous surfaces: 
a new look at Cassie's and Wenzel's laws 
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We consider a three dimensional liquid drop sitting on a rough and chemically heterogeneous 
substrate. Using a novel minimization technique on the free energy of this system, a generalized 
Young's equation for the contact angle is found. In certain limits, the Cassie and Wenzel laws, 
and a new equivalent rule, applicable in general, are derived. We also propose an equation in the 
same spirit as these results but valid on a more 'microscopic' level. Throughout we work under the 
presence of gravity and keep account of line tension terms. 
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I. INTRODUCTION 

Recently, there has been an upsurge of interest in 
wetting phenomena, particularly in the experimental 
arena with many dramatic new techniques and results 
available. The impact of these on theory has been 
twofold, both allowing older work to be thoroughly tested 
and also posing fresh challenges for those researching in 
the field today (for current reviews see B). Presently, 
most theoretical approaches are trying to go beyond 
wetting on simple, flat, homogeneous substrates and 
are exploring the effect a structured surface can have 
on the wetting phase. This heterogeneity can be both 
geometrical and chemical. However, there have existed 
for quite some time a few empirical laws describing the 
contact angle of a drop sitting on such a heterogeneous 
surface. Wenzel introduced an 'average' contact angle 
on a rough, chemically homogeneous substrate which 
is expressed in terms of the contact angle on a planar 
one. Likewise, for smooth but chemically heterogeneous 
surfaces the Cassie equation || is widely used, which 
defines an 'average' contact angle by the weighted mean 
of the angles that the drop would take on pure substrates. 
The associated wall tensions have been recently studied 
both for the geometrically rough [QJs) and for the 
chemically heterogeneous case || . Much of the literature 
has also been concerned with the modification of these 
laws when line tension effects are included (see, for 
example, [@-|9)). 

In this paper we investigate the statistical mechanical 
foundations of such empirical approaches. Young's 
equation JIo| is surely the bedrock of all wetting theory 
and from its most general form, after making an 
important number of assumptions, we are able to provide 
a systematic derivation of both the Cassie and Wenzel 
laws. We choose to describe the chemical heterogeneities 
within the surface in terms of interfacial and contact line 
tensions which are position dependent. We consider a 
solid substrate composed of several different atomic (or 
molecular) species and first coarse-grain up to a certain 



length scale, the small scale cut-off, in order to define 
appropriate composition variables. For the simplest case 
of a binary system, only one such composition variable is 
needed and can be defined, for example, as the relative 
area fraction of one of the two species in the surface 
layer. In this way, we arrive at composition variables, 
X(y) say. These, in general, depend on the coordinate 
y = (2/1,2/2) of the two dimensional surface and thus 
reflect the chemical heterogeneities which are present on 
length scales large compared to the atomic scale. 

A chemically homogeneous surface is characterised by 
position independent composition variables X(y) = X 
and one may define the different interfacial tensions a 
and the line tension A in the usual way JfjJ. The values 
of these tensions will, of course, depend on the values 
of the composition variables: a — a(X) and A = X(X). 
Thus, in the heterogeneous case with X = X(y), one may 
allude to a small gradient expansion and assume that 
the local tensions are given by a = a(X(y)) and A = 
A(X(y)), i.e. they are essentially determined by the local 
surface composition. In general, the anisotropy of the 
solid substrate will lead to a line tension which depends 
on the orientation of the contact line. In the following, 
this anisotropy will be ignored and the surface of the solid 
will be treated as a structureless wall. 

For imprinted surfaces, one has surface domains which 
are large compared to the small scale cut-off [ jL2| . 
An example is provided by domains obtained from 
microcontact printing which have typical sizes in the 
micrometre range. In this case, the various tensions 
are constant within the domains but vary across their 
boundaries. Droplets on such domains may exhibit 
contact angles which do not satisfy Young's equation 
in the limiting case where the boundary width is small 
compared to the domain size. However, in the absence 
of such an extreme separation of length scales, one has a 
position dependent Young equation as discussed in |T^] 
(though here contributions from the line tension have 
been ignored). 

The outline of our paper is as follows. In Sec. 
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H we specify the free energy of a drop adsorbed 
on a heterogeneous substrate. By applying a novel 
minimization procedure to this free energy, detailed in 
Sec. Q], we are able to find the most general form of 
Young's equation. A statistical mechanical interpretation 
of Cassie's and Wenzel's laws is then discussed in Sec. VI 
and a number of strong assumptions highlighted which 
are implicity adopted by them. In the following sections, 
using the generalized Young equation, we are then able to 
go on and provide a systematic derivation of both these 
rules where we keep account of line tension and gravity. 
For substrates which have both geometrical and chemical 
structure a new composite equation is suggested in Sec. 
IX . An alternative approach is adopted in Sec. [x] in which 
a local version of Cassie's law is proposed which needs 
no suppositions about the drop or heterogeneity of the 
surface. Finally, we make some short conclusions. 



II. WETTING OF A HETEROGENEOUS 
SURFACE 

We start by considering the free energy of a three 
dimensional drop of non-volatile [3 phase within a bulk 
a phase sitting on a rough, chemically heterogeneous 
substrate or wall. Let the area covered by the drop be 
r and let the edge of this area, i.e. the position of the 
contact line, be dT (see Fig. |l|). 




FIG. 1. Schematic of the three phase contact line of a 
droplet sitting on a rough, heterogeneous substrate. The 
contact area is denoted by T and the position of the contact 
line by dF. The upper right-hand quadrant denned by 
< y\ < Y is described by yi = 7(1/1) for some function 
7(2). 



The a(3 interface has a surface tension a and is at 
height h(y) above some reference plane, with Cartesian 
coordinates {y 1,1/2) = y as in Fig. jj The line tension of 
the contact or triple phase line, dT, is denoted A, while 
the wall-a and wall-/3 interfacial free energies are a wa 
and <j w p, respectively. If the substrate surface has a 
configuration Z{y) above the reference plane then the 
free energy is 



F[h, Z] = J dyl ay/1 + (Vh) 2 - \a wa (y) - a wfj {y) 



Vl + (VZ) 2 + Q(/i, Z) 



ds A(s) 



(2.1) 



dT{Z) 



Here s is the arc length of the contact line at dT on the 
substrate surface z = Z(y) and V is the two-dimensional 

gradient operator, V = (gf^af^)- We choose the 
location of the y-plane to be such that 



J dyZ(y) = 



(2.2) 



The function Q(h, Z) in ( |2.l|) can be decomposed into 
Q(h,Z) = Ap(h- Z) + ^Apg(h 2 - Z 2 ) (2.3) 

where Ap and Ap are the pressure and density differences 
between the (3 and a phases, respectively. The gravita- 
tional acceleration is denoted g. Intermolecular forces, 
such as van der Waals interactions, enter only implicitly 
through the various tensions |±3| . 



III. MINIMIZATION PROCEDURE 

The equilibrium configuration of the drop will be 
given by that location h(y) of the liquid surface or 
af3 interface and t hat contact line configuration dT(Z) 
which minimizes ( [2.l| ). To carry out this functional 
minimization it proves useful to take advantage of a 
version of Gauss's theorem in two dimensions. For future 
reference we state this here 

/ dy(V/).g = f figidyi - gidy 2 ) - / rfy/V.g 
Jr Jdv Jr 

(3.1) 

where the line integral is evaluated in a clockwise 
direction as indicated in Fig. |l|. 

Rather than functionally minimize with respect to h(y) 
with the variable boundary condition 



My) = z{y) 



for y G dT 



(3.2) 
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we choose instead a different approach. Initially dT 
is fixed, the equilibrium liquid surface calculated for 
this particular configuration and then the free energy 
minimized again with respect to the location of the 
contact li ne d T. 

Using ( |3.l|) minimization of ( [2.l|) with respect to h(y) 
is straightforward and we find the Laplace equation 



crV. 



V/i 



dQ 

+ -g{h,Z)=Q (3.3) 



with condition (3^). Note that this condition implies 
that the equilibrium h{y) will be a functional of dT. If 
we ignore gravity and the effective interface potential 
V, ( |3.3D implies that the surface has constant mean 
curvature Ap/2a. 

To consider variation with respect to the contact line 
we can without loss of generality specialize to the upper 
right-hand side quadrant of Fig. |l|. We assume that y on 
dT can be written as y = (2/1,7(2/1)) for some function 
7(1/1) and for < 2/1 < Y, where Y is defined by 7(F) = 
(see Fig. Q). Thus, the function 7(2/1) describes the shape 
of the contact line when projected onto the y-plane. 

For this region the free energy functional (2.1) can be 
written as (in the Monge representation) 



F[h,Z] 



dyi 



7(2/1) 



dy 2 { (7^1 + (V/i) 5 



(y) U/l + (VZ) 2 + Q(h; Z) 



dyi A (2/1, 7(2/1)) 



1 + 7W 



dy 



■Z(yi ,7(2/1)) 



(3.4) 



with the prime denoting differentiation with respect to 
the argument. The remainder of this part of the paper 
becomes by necessity quite mathematically intense and a 
casual reader may prefer to ski p im mediately to Sec. IV. 
Functional minimization of (|3.4|) leads to 



SF 

Sj(x) 



h 



h 



(3.5) 



where 

h = { o\J\ + (V/1) 2 - cr wa (y) - a wf 3(y) 



x^l+lYZ)- +QIZ.Z) 
S 



(3.6) 



Y 



dyi A (2/1, 7(2/1)) 



y=(x,~f(x)) 

l+7'(2/i) 2 



+ (^(2/1,7(2/0) 



(3.7) 



and 



h = 



07 



, , dQSh 



(3.8) 



The term ^3 arises from the implicit dependence of h(y) 
on "f(x) and can be derived using the chain rule for 
functional derivatives. 

We shall concentrate on simplifying the expression I 2 
first. This becomes 



h = A 



(l' Z yi ~ Z y 2 )( Z yiyi + 2 l' Z yiy 2 + 7 2 Z y 2 yi) 



- 7 "(1 + (VZ) 2 ) 



+7 2 (1 + C) 



2-f'Z yi Zy 



KA 1 + Z Vl Z V2 + Z y x ) 



-i\ yi {l + Z yi Z y Ji + Zl) 



+2 1 'Z yi Z y2+1 ' 2 (l + Zl) 



where the notation f x = j£ is used and all functions 



(3.9) 



are evaluated at y = (a;, 7(0;)). Fortunately, ( p.9[ ) can be 
re-written as 



J 2 = (VA.m - AC 9 )v/l + (V^) 2 



y=(x,7(x)) 



(3.10) 



where m = rii(y) is the unit vector orthogonal to both 
the normal to the surface z — Z(y) and to the tangent 
vector to the curve dT. Here we use the convention 
that the z-component of VA is zero. Furthermore, C g = 
C g (y) denotes the geodesic curvature of dT at y, i.e. the 
component of the 'acceleration' vector of the curve in the 
direction rh Hk 

Extending (|3.8| ) to the whole drop, one can use Gauss's 
theorem and write 



07 



/!= // y i V + (vM 2 + ^ 



Sh I hy^dyx - h yi dy 2 \ 
or h \ y/l + (Vhf~ ) 

dy { -crV. ; 

Sh I h y2 dyi ~ h yi dy 2 \ 



or 



Si \ ^l + (V/i) 2 




by virtue of the Laplace equation (3.3). Once again we 
now specialize to the up per right-hand quadrant. Here 
the boundary condition (ET^) can be expressed as 



h(yi, 7(2/1)) = -£(2/1,7(2/1)) 



(3.12) 



Functionally differentiating (3.12) with respect to 7(2;) 
leads to 
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dy 2 



x) + 



Sh 

Sj(x) 



OV2 



x) 



while differentiation with respect to y\ gives 

dh dh . dZ dZ . 

oyi oy 2 ayi dy 2 



(3.13) 



(3.14) 



As we are on dT, 7' (3/1) 



and so both (pU3|) and 



(3.14) allow simpl ificat ion of (3.11). Using (3.13) to 

Sh 

S~f 



eliminate 1^ and ( 3.14 ) to express everything in terms 



of yi , the integral can be evaluated and ( p.ll|) becomes 
Vh.V(h- Z) 



h = -o- 



(3.15) 



y=0,7(z)) 



IV. THE GENERALIZED YOUNG EQUATION 



We are now in a position to interpret ( 3.5|). At 



equilibrium ^ must vanish and so using (3.10) and 



( jyJl ), © implies 



-Vl + fVZ) 2 a wa (y) - a wlS (y) - VA.m(y) 



+A(y)C fl (y) 



(4.1) 



for y 6 dT. Recall that C g (y) is the geodesic curvature 
of the contact line at the point z — Z{y) on the substrate 
surface and m(y) is the unit vector perpendicular to the 
surface normal and to the vector tangential to the triple 
line at y. 

It follows from ( pT3| ) that Q(Z, Z) — 0, which we use 
to write (4.1) as 



1 + Vh.VZ 



o wa {y) - <J w p(y) - VA.m(y) + A(y)C g (y) 



(4.2) 



In fact, if we define the contact angle 6(y) to be the angle 
between the normal (and hence the tangent) vectors to 
the surfaces z = h(y) and z = Z(y) at y, then it is simple 
to show that 



cos[0(y)] 



1+VLVZ 



#+PF)(r+(W) 



(4.3) 



Hence, (4.2) is the generalized Young equation for a 
sessile drop sitting on a rough, chemically heterogeneous 
surface z = Z(y), 

cr wa (y) = CTioja (y) +ocos[0(y)] + VA.m(y) - A(y)C 9 (y) 

(4.4) 



where y is at the three phase contact line. This is the 
main result of our paper. 

Particular cases of the generalized Young equation 



(4.4) have been obtained previously. If the line tension 
terms are ignored, one has 



v wa {y) = o- w/3 (y) +crcos [6>(y)] 



(4.5) 



as derived by Lenz and Lipowsky 12|. 

The correction terms arising from the line tension have 
also been obtained previously for some special situations. 
First of all, their effect has been determined for planar 
and homogeneous substrates by Boruvka and Neumann 
p5| (and see also fTB]). This case will be discussed in the 



following section. In addition, Rusanov |17| has initiated 
study of an axially symmetric geometry for which the 
rotation axis is perpendicular to the y-plane. In this 
case, the shape Z of th e substra te surface depends only 
on the distance p = \Jy\-Vy\ from the rotation axis. 
Therefore, the wall tensions a w i = (J W i(p) and the line 
tension A = A(p) are taken to depend only on p and the 
problem becomes one dimensional. The contact line of 
the droplet term forms a circle of radius p = R, the unit 
vector rh has the p component m p = 1/yl + Z'(R) 2 , 
and the geodesic curvature is given by 



C g (R) = - 



1 



Ry/1 + Z'{Rf 



(4.6) 



If these expressions are inserted into the general 
equation (4.4), one obtains 



(Twa(R) = <J w p{R) + acos[9(R)} + 
x cos[4>{R)] 
where the slope angle <j> satisfies 



X(R)+X'(R) 



cos[c/)(R)} = 



1 



^l + Z'(Rf 



(4.7) 



(4.8) 



In order to understand the geometric meaning of this 
angle, consider the contour z = Z(p) of the substrate 
surface within the (p, z)-plane, and construct the straight 
line which is tangential to Z(p) at p — R. The slope angle 
4>(R) is the angle between this tangential line and the p 
axis. The special form ( |4.7| ) of the generalized Young 
equation (14.4! ) is equivalent to the equation derived by 
Rusanov |l7f| . 

Real substrates which are heterogeneous do not ex hib it 
the axial symmetry, assumed in the derivation of (|4.7|) , 
and the shape of the contact line will not be circular. In 
order to incorporate the deviations of the line shape from 
a circle in a qualita tive manner, Rusanov has proposed 
a generalization of (4.7) which involves a heuristic line 
roughness factor jl8] . 

Another effectively one dimensional geometry has been 
investigated by Marmur 
'cylindrical' interfaces whic 



|2fJ. He has considered 
1 depend only on one surface 
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coordinate, say y\ , and are translationally invariant 
in the indirection. For this case, the contact line 
is perfectly straight and lies at y\ = ±Y, for some 
constant Y say, implying that the geodesic curvature 
C g (±Y) = 0. The unit vector rh has a y\ component 
m yi = 1/ >/l + Z'{Y) 2 at 2/1 = y and so ( jOj ) becomes 

= a-^CT) + <7COs[0(Y)] + A'(Y)cos[^(Y)] (4.9) 

as stated in [E0| . 



V. PLANAR AND CHEMICALLY 
HOMOGENEOUS SUBSTRATES 



To show that (4.4) includes the more familiar form of 
Young's equation is quite straightforward. We specialize 
to a smooth, homogeneous substrate, that is 



A(y) = A 



0"uj/3(y) = °~w0 

z(y) = o 



(5.1) 



For planar curves the geo desic curvature simply becomes 
the curvature, looking at (|3.9|) the first term tends to 



AC= -A- 



7 "(x) 



(5.2) 



as Z(y) tends to a constant, while the vector m(y) in the 
s am e limit is the normal n(y) to the curve. Consequently, 
(4.4) is now 



0~w/3 + & COS 9 — AC 



(5.3) 



This is Young's equation [[To| with the Boruvka and 
Neumann line tension term. 

If corrections due to gravity and long-range inter- 
molecular forces are ignored, the solution of the Laplace 
equation ( |3 .3| ) is a spherical cap with radius of curvature 
i?7r = 2er/Ap. It is then easy to see that the curvature of 
the contact line itself satisfies 



C 



1 



sin 9 
Ap 
2a sin 9 



(5.4) 
(5.5) 



When this relation is inserted into (5.3), the contact angle 
is found to depend on the two dimensionless parameters 
{&wa ~ o~ w p)/cr and ApX/a 2 . Note that, in the pressure 
ensemble considered here, the curvature C of the contact 
line is fixed for a smooth, homogeneous substrate. 

Similarly, if one specifies a fixed volume V of j3 
phase then again for a smooth, homogeneous substrate 
a relationship between C, 9 and V can be determined. 
Ignoring gravity and long-range forces, the droplet must 
take the form of a spherical cap. C onse quently, it is easy 
to determine its volume and using (5.4) relate this to the 
contact line curvature 



G 



7r(2 - 3cos6» + cos 3 



3ysin d 



1/3 



(5.6) 



The contact angle once more depends on (o~ wa — cr w ^)ja 
but now also on the new dimensionless quantity X/aV 1 ^ 3 . 
In principle, the line tension A may be determined via 



( |5.3| ) if one knows the interfacial tensions a wa , o w p and 
a and if one measures the geometric parameters and 
C. The contact line curvature C depends on both in 
the Ap- and in the V- ensembles. For the latter, it also 
depends on the droplet volume V, see (5 



Finally, thermally excited fluctuations away from the 
minimal free energy would in general change the contact 
line contour and so give additional corrections to 0. Such 
effects are not included in this paper. 



VI. SPATIAL AVERAGING 



Using the generalized Young equation (4.4) it is 
possible to provide a systematic derivation of the 
Cassie law (valid for flat but chemically heterogeneous 
substrates) and the Wenzel rule (applicable for purely 
rough surfaces), 



cos 6> e ff = /lfi cos 0i 



cos c s — r cos 



Cassie (6.1) 
Wenzel (6.2) 



Here 9 e g is an effective or average contact angle assumed 
on the heterogeneous substrate, 9i is the angle taken 
on a simple planar surface composed entirely of surface 
component or chemical species i, fi is the fraction by area 
of the surface made up of i and r is the ratio of true to 
planar surface area. 



The starting point for a derivation of equations (6.1) 
and (6.2) should be a definition of the effective or average 
angle 9 Q s . However, in order to perform any average one 
must first establish choices for the following criteria: 



(i) what statistical mechanical ensemble should be 
used and what are the different states in this 
ensemble? 

(ii) what are the a priori probabilities or statistical 
weights assigned to these different states? Often, 
one makes an implicit assumption about equal a 
priori probabilities but this choice is not unique. 

(iii) which quantities does one wish to average? 

To answer question (i) we can choose to prescribe 
either the volume of the f3 drop or else the pressure 
difference, Ap, across the a(3 interface. It is more 
convenient to use Ap as the basic variable and so we 
opt to work in the corresponding pressure ensemble. 
Consequently, by ignoring gravitational and long-range 
intermolecular force corrections, we have droplet surfaces 
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of constant mean curvature. For flat and homogeneous 
substr ates , the latter are given by spherical caps and, 
from ( |5.5| ), have a fixed contact line curvature C. As 
discussed before, this is no longer t rue in the volume 
ensemble since C depends on V, see (5.6). Therefore, if 
one wants to deduce the magnitude of the line tension A 
from observations of the droplet shape, it is important to 
specify the appropriate ensemble one is working in. 

Having selected an ensemble we next need to decide 
what the different states are in this ensemble. Three 
different possibilities present themselves: 

(a) One could investigate a drop at a certain position 
on the surface and move along its contact line. 
Aside from the flat and homogeneous case, the 
contact angle will vary and thus, one could study 
the average 9 of an individual droplet. If a flat 
and homogeneous substrate is consi dere d then the 
contact angle is constant satisfying (5.3). 



(b) Different positions for the drop could be considered 
and the different states would then be determined 
by these positions. For the flat and homogeneous 
case this is the same as (a). However, if one places 
a drop at different positions on a heterogeneous 
substrate, the contact line contour and the contact 
angles at that contour will, in general, change if we 
require the state of the drop to be a local minimum 
of its free energy. The task of determining these 
states and then assigning a priori probabilities to 
them seems to be prohibitively difficult. 

(c) Instead of considering the states of the drop one 
could instead opt for those of the contact line. 
To start, a certain position y on the surface is 
selected and an average over all triple line contours 
which pass through this point and correspond to 
the states of the droplet which are local minima of 
the free energy is carried o ut. To be local minima 
these states must satisfy (^3), and so (if gravity 
and long-range intermolecular interactions can be 
ignored) be of constant mean cur vature, as well as 
obeying the Young equation (4.4). A final further 
average over all positions y is then taken. 

In this section of the paper we choose option (c) 
above and proceed to average over different contact lines. 
To answer (ii) a Boltz man n weight is chosen for each 
contour, calculated via (2.1), and again we point out that 
only those contours which are local minima of the free 
energy a re to be considered. 

From (4.4) one can see that the variation of the local 
contact angle 8(y) is mainly governed by the change 
in the interfacial free energies a wa and a w p (the line 
tension is typically » 10 _9 Jm _1 ||). However, these 
quantities determine the cosine of the contact angle and 
consequently summing over their position dependence 
should give an effective cos£? e ff- Therefore, in response 
to (iii), we choose to average cos[0(y)] rather than 6(y) 



itself. Arguing in terms of surface tensions, it is the 
component of the a(3 tension that is in the surface that 
is the one chosen to be averaged. 

The Young equation (4.4) implies that the contact 
angle is dependent on position y and, as mentioned 
before, the local configuration taken by the triple line 
contour dT at y, that is 



e = 6(y;dT) 



(6.3) 



To proceed, we define {dr(y)} to be the set of all 
contours of the drop which pass through y and are 
local minima of (2.1). Then option (c) can be written 
mathematically as 



COS Or. 



IdyE {d r {y n^(e(y;dT))e- F [ h ^ z V T 



'{ar(y)> 
JdyE 



{dr(y)}' 



-F[h[dF],z]/T 



(6.4) 



where we use a Boltzmann weight function with units 
such that fcg = 1 and write explicitly the functional 
dependence of the equilibrium h = h[dT] on the contour 
configuration. 

To make further progress several strong assumptions 
are needed. Firstly, the dependence of the contact angle 
on the shape of the droplet is ignored, i.e. 



cos[0(y;ar)] «cos[0(y)] 



(6.5) 



Secondly, the heterogeneities are taken to be such that 
the drop is not strongly confined to or repelled from any 
particular region. All non-homogeneities are 'uniform' in 
this aspect. The probability that the contact line of the 
droplet passes through y is given by 



E 



{9r(y)}' 



'[fc[ar],z]/a 



JdyE< 



-F[h[dr],z] IT 



(6.6) 



<{ar( y )}' 

and we assume that this is approximately independent of 

y, 



E 



{£>r(y)} ' 



-F h[dr] ,Z /T 



1 



1 



{9r(y)} ' 



-F[h[dr],z]/T J dy A.* 



(6.7) 



where A n represents the projected area of the substrate 
surface onto the y-plane. No posit ion is dramatically 
favoured over any other. Equation (6.7) clearly drasti- 
cally r edu ces the domain of validity of any results derived 

from (6.4). Notice, th at for a flat and homogeneous 

substrate (|6.5| ) and (6/7) are satisfied identically. 

Consequently, (3.4) reduces and the effective angle now 
satisfies 



cos 6> off 



Jdycos[e(y)} 

An 



(6.8) 



which we use as a working definition. 
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VII. PLANAR BUT CHEMICALLY 
HETEROGENEOUS SURFACES 

We first specialise to smo oth, chemically hete rogeneous 
substrates. In the light of (6.5), equation (4.4) is simply 



<7wa(y) = Otu/j(y) +CTCOs[0(y)] 



(7.1) 



where all the line tensions terms are ignored. When this 
is integrated over the surface we are led to (recall (|6.8|)) 



A^ct cos 6> off = / dy[a wa (y) - o- w0 (y)] 



(7.2) 



where A v — J dy. If our substrate has a fraction (by 
area) /, of substance i, then 



dy [cr wa (y) - cr w p(y)] = Ar / , fifota ~ <?ip) (7-3) 



where o~i a is the wall-a interfacial free energy for 
a substrate c hem ically coated with material i, etc. 
Consequently, mTa) is 



COS0 eff = ^ fi( a ia ~ °ip) I <* 



(7.4) 



which can then be interpreted in terms o f mo re physical 
parameters using the Young equation (5.3). At this 
point, we would like to emphasize that (7.4) already 
contains, at the level of the approximations used (see 
(6.5) and (6.7)), all the physics associated with the 



heterogeneous substrate which is now completely char- 
acterised by the numbers /j. 

Therefore, if, in the spirit of Cassie, we use 



a lfj = a cos 9 i - A, Q 



to re-express (7.4) in the traditional form 



cos tf cff 



COS 6; 



(7.5) 



(7.6) 



then the values of Bi and Ci will, in general, have no 
relation to the original experimental set-up. The Young 
equation, (7.5), is obeyed by all droplets adsorbed on a 
chemically homogeneous surface i. Therefore, provided 
6i and Ci are measured for the same drop, it does not 
matter whether the pressure or volume ensemble is used, 
or what the pressure difference or volume of (3 actually 
is. If one chooses to make the measurements on a 
heterogeneous substrate then it is vital that a point is 
selected for which the chemical composition of the surface 
can be unambiguously determined as i, and that both 
the contact angle and curvature of the contact line are 
measured precisely at this point. 



We feel that this step, from ( [7.4] ) to (|7.q), has not been 
highlighted in the literature, and, while it is certainly 
possible to determine Ci and 6i from a particular 



point of a complicated drop configuration present on a 
heterogeneous surfac e, t here is no advantage in doing so. 
We emphasize that (7.4) is the fundamental equation in 
which no ambiguities arise. 

Note that the line tension terms in ( |7.6D are identical 
to those first included by Drelich and Miller Q by the 
use of thermodynamic arguments. 



VIII. CHEMICALLY HOMOGENEOUS BUT 
ROUGH SUBSTRATES 

Wenzel's rule is taken to be valid for rough, chemically 
homogeneous substrates. However, from (4.4), when line 
tension terms are ignored, the contact angle taken on a 
rough surface (defined as the angle between the normal 
vectors to the surface and to the drop) is identical to 
that on a planar substrate. Consequently, we believe 
that (6.2) in fact refers to the average of a local planar 
contact angle ^(y), defined as 



cos[6' 7r (y)] 



1 



(8.1) 



year 



Thus, 9-rr(y) is the angle between the tangent to the 
drop's surface and a local horizontal plane given by 
r=(y',Z(y)) for all y'. 

Such an approach has a long history with, we believe, 
Shuttleworth and Bailey ]2lJ first defining the observed 
contact angle to be the sum of the actual (or in their 
terminology intrinsic) contact angle and the slope angle, 
(see Fig. |) 



My) = o(y) + <t>(y) 

Here the slope angle is defined by 

1 



cos[c/>(y)] 



V 7 ! + NZ(yW 



(8.2) 



•3) 




FIG. 2. A cross section through a drop sitting on a one 
dimensional r ough surface z — Z(y). Th e tru e contact angle 
0, defined in (4.3), the local angle, 9-n, (8.1), and the slope 
angle cf>, (B.5), are all clearly marked. 
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For a substrate surface shape z = Z(x), which 

is dependent on only one surface coordina te x, it 

is straightforward to show that ( |S.l| ) and ([T^) are 
identical. For example, if the substrate surface is axially 
symmetric and its height is des cribed by z — Z(r), we 
obtain the slope angle </>, (O), which has the simple 
geometric interpretation explained in Sec . |rv| . In a three 
dimensional space the decomposition (^^) is not valid in 
general as there is no guarantee that <fi and 9 lie in the 
same plane. Our choice, 9^ or the slope angle of the a(3 
interfacial surface, is more general and does not suffer 
from any such problems. 

Proceeding wi th t he definition ( |6.8| ), valid for 9(y) = 
On{y)> we write (4.4) as 



tr cos[0,r(y)] + a 



Vh.VZ 



(8.4) 



where (3.5) is again called upon to ignore line tension 
effects. The surface area of the rough substrate, A, is 



a = J dy^/i + (yzj 



(8.5) 



and therefore from ( |3.l| ) and (p.3|), the integral of ( |8.4| ) 
becomes 

f dQ 

AkV cos 9 cS - a I dyZ(y) — (Z, Z) = A(a wa - a wl3 ) 

(8.6) 



IX. ROUGH AND CHEMICALLY 
HETEROGENEOUS SURFACES 



Continuing to use the definitions (6.8) and (3.1) for the 



effective contact angle, it is not too difficult to see that 
an additional law is possible for a sessile drop on a rough 
and che mically heterogeneous substrate. Equation (|4.4|), 
via (6.5), 



<rcos[0 7r (y)] +(T 



Vh.VZ 

\A+( V ' 1 ) 2 



(y) - Vwp{y)]\jt + (VZ) 2 

can be integrated over y. Using (|3.3|) one finds 



(9.1) 



f dQ 
Ana cos 6 e g = ex / dyZ{y) — {Z,Z) 



dy 



(y) - °wf3 (y) 



Vi + (v^) 2 



and so arrives at 
cos 6> cff = ^ '< 



ApgZ 2 



(9.2) 



(9.3) 



where is the ratio of the non-planar surface area 
covered with material i to the total planar area. 



X. AN ALTERNATIVE PRESCRIPTION 



cos0 eff - -j- Jdy (A P Z(y) + ApgZ 2 (y)) 



•7) 



using ( p.3j ) and (p.3[), where 9 and C can again be 
measured for any drop located on a planar substrate. 
Notice that Z(y) — > as y — > oo has been assumed. 
Denoting r = to be the ratio of non-planar to planar 
surface areas and utilizing the definition ( [2.2j ) we have 



cos 9 cS 



cos6» ) +ApgZ' 2 

a 



where Z 2 is the mean square height of the surface 
^_ fdyZ 2 (y) 



(8.9) 



Equation ( p.8| ) is Wenzel's rule incorporating line 
tension effects and gravity for the first time. Here there 
is some disagreement between our line tension correction 
and that proposed by Drelich It is also interes ting 



to note that the simple form of Wenzel given by (6.2) 
was postulated a long time ago to become invalid under 
a gravitational field p2[ . Finally, we again wish to 
emphasize that (3.6) is independent of the properties of 
the drop being considered. 



It is somewhat unsatisfying that the traditional forms 
of the Cassie and Wenzel equations are independent of 
the shape of the drop and assume a ' unif orm' distribution 
of the heterogeneities as given by (6/7). In this section 
we derive a new relation for the effective contact angle 
9 e g on a chemically heterogeneous substrate, which while 
not being quite as aesthetically pleasing, does not have 
these two drawbacks. Here we switch from the choice 
(c), described in Sec. VI, of the states in our statistical 



mechanical ensemble to that of (a). 

Let a drop of (3 phase on a flat, chemically heteroge- 
neous substrate equilibriate and take up the shape of the 
optimum contour, i.e. that which is a local minimum of 
the free energy. We then define the effective contact angle 
to be simply the average contact angle taken around this 
particular configuration of the three phase contact line. 
Writing s for the arc length on dT and the local contact 
angle then as 9(y) — 9(s) we have 



cos 9 ctf = 



<f gr ds cos [9(s) 
§ dr ds 



(10.1) 



Integrating Young's equation in the form of (7.1) around 
the optimum contour gives 
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a cos 9 c ff (p ds — 



ds 



or 



(s) - <J w p{s) 



+ 



ds 



or 



X(s)C(s) - VA.n(s) (10.2) 



Now looking again at the upper right-hand quadrant of 
Fig. [l], the normal to the contact line is given by 



n(y) 

and so 

V.n(y) = - 

implying 
Hence, we have 



Vl + 7'G/i) 2 



1 



(10.3) 




V.n(s) = -C(s) 



-c(yiMyi)) 

-C(y) (10.4) 
(10.5) 



cos 9 c g 



U {p ia — (Jiff) 



§ dr dsV.(\(s)h(s)) 



far 



ds 



(10.6) 



cos 9 cS = 



cos 6*, 



$, r dsV.(A(s)n( s )) 



far 



rf.s 



(10.7) 



where £i is the fraction of the total perimeter of the drop 
which is on surface composed of material i and A^ is the 
line tension for this surface. The angle 0i and curvature 
Ci are defined as before and can be measured for any 
drop on any flat surface i providing both measu rements 
are taken at the same point. Equation (10.7) can be 
written quite succinctly as 



cos( 



eff = cos ( 



(10.8) 



if line tension contributions are negligible. 

Fr om a theoretical point of view (10.7) is preferable to 
(ph as 



(i) only one (the most likely) contour is considered 

(ii) it is a local equation and it is only the surface in 
the immediate neighbourhood of the contact line 
that is need ed. It would be obviously foolhardy 
to use (7.6) for a surface half of which is purely 
hydrophobic, for example, and the rest consisting 
of alternating hydrophobic and hydrophilic strips if 
the drop sits entirely in either of these two regions. 
Equation ( |10.7 ) requires no assumptions about the 
nature of the heterogenieties. 



(iii) there is explicit dependence on the droplet shape 
due to the last integral term (remember that n is 
the normal to the contact line) 

Experimentally, during the last decade a variety of 
techniques have been developed in order to measure sur- 
face topographies on small scales. One pertinent example 
is provided by atomic or scanning force microscopy, which 
makes it possible to measure the precise shape of small 
droplets on nanometre scales and so, the exact location of 
the contact line, see e.g. [p3[ . Similar techniques can also 
be used in order to determine the chemical composition 
of the substrate surface. Consequently, one may obtain 
an estimate for the variation of the line tension A along 
the contact line dr. T herefore, it seems probable that the 
line integral in (10.7) can be determined experimentally 
from a careful analysi s of scanning force micrographs. 
Hence, we can expect (10.7) to be of some practical use. 



XI. CONCLUSION 

In this paper, by introducing a novel minimization 
scheme of the free energy, we have found the most general 
form of the Young equation. Equation (4.4) is valid for 
substrates both chemically and geometrically inhomo- 
geneous, in the presence of long-range interactions and 
under the influence of gravity. Taking advantage of this 
new result, we have looked again at the phenomenological 
laws of Cassie and Wenzel and examined their statistical 
mechanical foundation. 

To recover ( |6~T ) and (|6.2|), two strong assumptions are 
needed: (i) the dependence of the contact angle on the 
shape of the drop has to be ignored and, (ii) all locations 
of the droplet are taken to be equ ally likely. These two 
assumptions lead to the estimate (3.8) for the effective 
contact angle. 

A new insight, that our approach has brought to 
light, is the interpretation of the contact ang le 0i and 
triple line curvature Ci in the Cassie, ( |7.6[ ), and in 
the Wenzel, 



(1 



equations. Recall that here the 
index i distinguishes between the different compounds 
or chemical species in the substrate surface. When using 
these results it is vitally important to realize that di and 
Ci are present in the exact com bination described by the 
planar Young's equation, (5.2), i.e. for each substrate 
component i, a cos 9 i — XiC'i always appears, which is 
identically equa l to Oi a — cuq. Therefore, when we wish 
to understand (^^) and (|8.8| ) on a practical level, the 
effective contact angle is determined by the ratios fi 
or r and by the value of <7 ja — Oi@ for the surface i 
under inve stig ation. Due to the universality of Young's 
equation, (5.3), the actual 9i and Ci used are irrelevant, 
in the sense that they need bear no relation to the droplet 
shape observed in our current experiment. All that is 
required, is that they are measured at the same position 
in space for any drop of in a phase, providing it lies, 
at that point, on a surface domain composed of the 







substrate compound i. Young's equation then guarantees 
that the measured values of Oi and Cj, when used in 
the Cassie or Wenzel laws, will simply combine to give 
0m — cri/3 as require d. T o reiter ate , due to the particular 
way they occur in (7.6) and (8.8), 9i and Cj have, in 
general, no relation to the actual shape of the droplet on 
a heterogeneous substrate. 

Within this remit, the relation ( |9.3| ) is proposed which 
relates the effective contact angle on a substrate with 
both geometric and chemical heterogeneities to those on 
s impl e planar surfaces. We also find the new relation 
(10.7) for the effective contact angle on a chemically 
heterogeneous substrate, which is true generally and 
requires none of the above assumptions. Consequently, 
the influence of the surface on the shape of the wetting 
droplet is taken into account. As discussed in the 
previous secti on, th e line integrals which appear in the 
new relation ( |10.7| ) for the effective contact angle can 
be estimated from scanning force microscopy data of the 
droplet shape. Such an alternative approach does not 
recover a local form of Wenzel's rule, possibly because 
this rule does not involve an average of the true contact 
angle but rather of a local planar one. 
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